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Abstract 

Deuteron-deuteron elastic scattering and transfer reactions in the energy regime above four-nucleon breakup threshold are described 
by solving exact four-particle equations for transition operators. Several realistic nuclear interaction models are used, including the 
one with effective many-nucleon forces generated by the explicit A-isobar excitation; the Coulomb force between protons is taken 
into account as well. Differential cross sections, deuteron analyzing powers, outgoing nucleon polarization, and deuteron-to-neutron 
polarization transfer coefficients are calculated at 10 MeV deuteron energy. Overall good agreement with the experimental data is 
found. The importance of breakup channels is demonstrated. 
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1. Introduction 

The pursuit for numerical solutions of the three- and 
four-nucleon scattering problems has been, since the early 
seventies of the last century, one of the most challenging 
endeavors in nuclear reaction theory, following the develop¬ 
ment of formal exact N-body equations using momentum 
or configuration space representations [1,2]. Progress ad¬ 
vanced slowly at first and limited to the use of separable rep¬ 
resentations of subsystem operators, but with the advent of 
powerful computational tools, both in terms of algorithms 
(spline interpolation. Fade summation, special integration 
meshes and weights, etc) and hardware advances, three- 
nucleon (3N) calculations with realistic nucleon-nucleon 
(NN) force models reached state-of-the-art status in the 
early 1990’s due to the effort of a number of independent 
groups [3-7]. Due to its higher dimensionality and multi¬ 
channel complexity, the four-nucleon (4N) scattering prob¬ 
lem took twenty years longer to reach the same status 
as the three-nucleon system except for the calculation of 
breakup amplitudes. These developments are mainly due 
to the works of the Pisa [8-11], Grenoble-Strasbourg [12- 
15], and Lisbon [16-19] groups. Because the first two groups 
use the coordinate-space representation, they were able to 
include, not only realistic two-body interactions, but also 
realistic three-body force models. Nevertheless they have 
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had a major difficulty in calculating multichannel reactions 
and going beyond breakup threshold, particularly when 
the Coulomb interaction is included between protons. The 
Lisbon group uses the momentum space Alt, Grassberger 
and Sandhas (AGS) equations for transition operators [2] 
that can be solved for multichannel reactions both below 
and above breakup and with the Coulomb force included. 
The only stumbling block has been the inclusion of irre¬ 
ducible three-body forces. As alternative the nuclear force 
model with explicit excitation of a nucleon to a A isobar 
was used. This coupling generates both effective three- and 
four-nucleon forces (3NF and 4NF) that have been success¬ 
fully included in 4N calculations by the Lisbon-Hannover 
collaboration [20]. The calculations using potentials derived 
from chiral effective field theory have been performed as 
well [16,17] but so far including only the NN part of the 
interaction. 

In the last 40 years progress in nuclear reaction theory 
has most often succeeded experimental developments to the 
point that, when calculations achieved a solid ground, the 
instrumentation that gave rise to the data was no longer 
in operation. Therefore inconsistencies between different 
data cannot anymore be resolved by repeating the exper¬ 
iments or developing new ones guided by the theoretical 
predictions. The 4N scattering problem has suffered from 
this much more than the 3N system for the reasons men¬ 
tioned above. Nevertheless, new 4N scattering calculations 
are worth pursuing because they lead the way to the solu- 





tion on complex multiparticle scattering problems, not just 
in nuclear physics but also in cold atom physics [21]. 

In this work we present first results for 4N reactions 
initiated by the collision of two deuterons (d) at energies 
above four-particle breakup threshold. In this energy do¬ 
main there are a few shallow resonances [22]; therefore one 
does not expect the same problems as encountered in n-^H 
and n-^He near threshold. However, in two aspects the the¬ 
oretical description of the d + d scattering is even more in¬ 
teresting and challenging. First, since deuterons are loosely 
bound and spatially large systems, the scattering of two 
deuterons looks like the collision of two identical halo nu¬ 
clei. The coupling to breakup channels in such a system is 
considerably stronger than in nucleon-trinucleon scatter¬ 
ing. This makes the d + d reactions computationally more 
difficult since open breakup channels lead to most compli¬ 
cated singularities in the kernel of scattering equations and, 
furthermore, a larger number of partial waves is needed. 
Second, the deuteron being a spin one particle also pro¬ 
vides the opportunity to calculate a number of tensor ob¬ 
servables, both md + d elastic as well as in ^H((i,p)^H and 
^H(ci, n)^He transfer reactions for which there is experi¬ 
mental data. The most abundant set of the experimental 
data ior d + d reactions exists at deuteron energy Ed = 
10 MeV where not only differential cross section and an¬ 
alyzing powers but also deuteron-to-neutron polarization 
transfer coefficients have been measured [23] establishing 
the n)^He reaction as an efhcient source for polarized 
neutrons. 

In Section 2 we explain the AGS equations we use and 
how to solve them. In Section 3 we show results for d -|- d 
elastic scattering, while in Section 4 results for ^H(d,p)^H 
and ^H(d, n)^He transfer reactions are presented. Conclu¬ 
sions are drawn in Section 5. 

2. Deuteron-deuteron scattering equations 

As in our previous studies of 4N scattering [16,24,19], we 
take advantage of the isospin symmetry and treat protons 
and neutrons as identical fermions. This enables the sym- 
metrization of the transition operators W/Jq- thereby reduc¬ 
ing the number of components that are distinct according 
to two-cluster partitions [16]. As usual, a = I labels the 
3-I-I partition (12,3)4, while a = 2 stands for the 2-1-2 par¬ 
tition (12) (34). To ensure the required full antisymmetry 
of the four fermion system, the employed basis states have 
to be antisymmetric under the exchange of two particles 
in subsystem (12), and also in the subsystem (34) for the 
2-1-2 partition. The 4N transition operators 14/32 describing 
deuteron-deuteron scattering obey the AGS integral equa¬ 
tions 

^12 = (GotGo) ^ — P34UiGotGol4i2 + U2GotGol422, 

(la) 

U 22 = (1 - P34)UiGotGoUi2- (lb) 

Here P 34 is the permutation operator of particles 3 and 4, 


t = V + vGot (2) 

is two particle transition operator for the pair (12) inter¬ 
acting via potential v, while Ui and U 2 are symmetrized 
transition operators for 3-1-1 and 2-1-2 subsystems, respec¬ 
tively [16]. The dependence of all transition operators on 
the available energy E arises through the free four-particle 
resolvent 

Go = (E + ie - Ho)-\ (3) 

with EIq being the free Hamiltonian. The finite imaginary 
part is is introduced in the complex energy method to avoid 
singularities in the kernel of AGS equations, but the limit 
e -bO is needed for physical amplitudes, that are given 
by the on-shell matrix elements of the transition operators 
U/da [16]. The £ —> -bO limit is obtained by the analytic 
continuation of the finite e results using the point method 
[25]. The analytic continuation, however, is only accurate 
when using sufficiently small e values at which the kernel 
of the AGS equations, although formally being nonsingu¬ 
lar, still shows a quasisingular behavior [18]. These qua¬ 
sisingularities reflect the presence of open p + ^H, n -b ^He, 
d -b d, d -b n -b p, and n -b n -b p -b p channels. Their treat¬ 
ment is taken over from Ref. [18] where a special method 
for numerical integration absorbing the quasisingular fac¬ 
tor into the integration weights was developed. This way 
the quasisingularities can be integrated accurately without 
increasing significantly the number of grid points. In the 
present calculations we obtain well converged results using 
£ between 1.2 and 3.0 MeV in 0.3 MeV steps and about 30 
grid points for the discretization of each momentum vari¬ 
able. Note that due to a larger weight of breakup channels a 
bit more £ values (about 6) are needed for a reliable £ —>■ -bO 
extrapolation as compared to previous calculations [18,26]. 

Although we explore the isospin symmetry, we also ac¬ 
count for the isospin violation effects due to the pp Goulomb 
repulsion and the hadronic charge dependence (GD) of the 
nuclear force. These effects cause the two-nucleon transi¬ 
tion matrix t to couple the states with different total isospin 
in both 3N and 4N systems. In 3N or 4N total isospin 
basis, the two-nucleon transition matrix t is given by lin¬ 
ear combinations of pp, np, and nn transition operators as 
described in Ref. [27]. For the pp pair beside the nuclear 
force also the screened Goulomb potential is added, en¬ 
abling rigorous inclusion of the Goulomb interaction in the 
deuteron-deuteron scattering via the method of screening 
and renormalization [28,29,24]. We obtain fully converged 
results calculating the Coulomb-distorted short-range part 
of the amplitudes with the screening radius i? = 13 fm. 
The direct unscreened Coulomb amplitude is present only 
in the elastic scattering; it is added after the renormaliza¬ 
tion of the short-range amplitude. The direct Coulomb am¬ 
plitude causes the d-\- d elastic differential cross section to 
diverge in the forward and backward direction but is absent 
for transfer reactions that are only distorted by Coulomb 
[24]. Other electromagnetic effects as the magnetic moment 
interaction are not explicitly included in the calculations, 
however, their short-range part is implicitly included in the 



employed NN potentials that are fitted to the NN data. 
They are not isolated in the present work but, based on 
previous studies [17,30] where they have been found to af¬ 
fect the vector analyzing powers only and to decrease with 
increasing energy, at 10 MeV one could expect only minor 
effects. 

We solve the AGS equations in the momentum-space 
partial-wave framework following the methodology devel¬ 
oped in Refs. [16,18]. In this framework the AGS equations 
constitute a large system of coupled integral equations in 
three continuous Jacobi momentum variables. The inte¬ 
grals are discretized using Gaussian quadratures with spe¬ 
cial (standard) weights for quasisingular (nonsingular) in¬ 
tegrands, leading to a huge system of linear algebraic equa¬ 
tions. Beside the number of grid points that can be kept 
moderate, the size of the resulting system depends also on 
the number of included angular momentum states. With 
respect to the number of partial waves needed for achieving 
convergence, the d + d reactions are more demanding than 
those initiated by the p-|-^Horn-|-^He collisions already 
calculated in Refs. [19,27]. In terms of the basis states de¬ 
fined in Refs. [31,27], when solving the AGS equations we 
include 4N partial waves with orbital angular momenta lx 
up to 6 and up to 7, total angular momenta of the 
2N subsystem jx^jy up to 6, total angular momentum of 
the 3N subsystem Jy up to and total 4N angular mo¬ 
mentum J up to 7. Once the AGS equations are solved, 
for the calculation of elastic observables it is sufficient to 
include only the initial and final d d states with G < 4. 
In contrast, transfer reactions ^H(d,p)^H and ^H(d, n)^He 
require G at least up to 6 in the channel states. 

3. Elastic scattering 

As in our previous calculations of nucleon-trinucleon 
scattering, we use several realistic nuclear force models, 
enabling us to study the sensitivity of the predictions to 
the dynamic input. Beside two purely nucleonic interaction 
models, the inside-nonlocal outside-Yukawa (INOY04) po¬ 
tential by Doleschall [32,13] and the CD Bonn potential 

[33] , we also use the two-baryon potential CD Bonn + A 

[34] that is the coupled-channel extension of CD Bonn, 
explicitly allowing virtual excitation of a nucleon to a A 
isobar and thereby yielding mutually consistent effective 
three- and four-nucleon forces. This model, however, is 
not fitted to the trinucleon binding energy (BE), yielding 
7.53 (8.28) MeV for ^He (^H) which is increased relative 
to the CD Bonn BE result of 7.26 (8.00) MeV. Only the 
INOY04 model, predicts the BE of ^He (^H) to be 7.73 
(8.49) MeV, nearly reproducing the experimental value of 
7.72 (8.48) MeV. Since the p -I- and n J- ^He threshold 
positions depend on the respective BE, some scattering 
observables are expected to correlate with the BE, thereby 
establishing INOY04 as a reference potential. Of course, 
the dependence of the observables on the used interaction 
model is in general much more complicated, but in partic- 



Fig. 1. Differential cross section ol d + d elastic scattering at 10.5 
and 12.3 MeV deuteron energy. Results calculated using INOY04 
potential are compared with experimental data from Refs. [35] (□), 
[36] (a), and [37] (x). 

ular cases simple correlations with nonmeasurable bound 
state properties such as the deuteron D-wave probability 
Pd may take place. For the INOY04 potential Pd = 3.60 
% while CD Bonn and CD Bonn -|- A have Pd = 4.85 %. 

In the present work we concentrate on the d + d scat¬ 
tering at deuteron energy Ed = 10 MeV where the most 
abundant set of experimental data exists, both for elastic 
and transfer reactions. The differential cross section da/dft 
for elastic d + d scattering is, however, an exception. We 
therefore show da/dfl as a function of the center of mass 
(c.m.) scattering angle 0c.m. in Fig. 1 at Ed = 10.5 and 12.3 
MeV. Solely the INOY04 potential is used for predictions 
but, based on the study at Ed — 10 MeV, the sensitivity 
of this observable to the force model is small. The angu¬ 
lar dependence of da/dVL is simple, with forward and back¬ 
ward peaks, where da/dU diverges due to the long-range 
Coulomb amplitude, and a minimum at 0c.m. = 90°. This 
shape remains almost constant while the absolute value of 
the differential cross section decreases with increasing en¬ 
ergy. Regarding the description of the experimental data, 
the picture is a bit contradictory. At Ed = 10.5 MeV there 
is a good agreement with the data from Ref. [36] while the 
data from Ref. [35] are slightly overpredicted, by 3% at 
0c.m. = 90°. At Ed = 12.3 MeV the data from Ref. [37] 
are well reproduced by the calculations at 0c.m. < 65° but 
slightly underpredicted around the minimum, by 4% at 
0c.m. = 90°. These findings suggest that more calculations 
over a wider energy range need to be performed and com¬ 
pared with the available data to determine the discrepan¬ 
cies between theory and experiment and find out possible 
inconsistencies between data sets. 

In Fig. 2 we present results for deuteron vector analyz¬ 
ing power iTii and tensor analyzing powers T 20 , T 21 , and 
T 22 ind + d elastic scattering at Ed = 10 MeV. The predic¬ 
tions are obtained using the potential models INOY04, CD 
Bonn -I- A, and CD Bonn. These spin observables are very 
small in their absolute value, of the order of 0.02. Due to 
the identity of the two deuterons the angular distributions 
of elastic observables in the c.m. frame are either symmet¬ 
ric {da/dVL, T 20 , T 22 ) or antisymmetric {iTn, T 21 ) with re¬ 
spect to 0c.m. = 90°. The overall description of the exper¬ 
imental data is good. The data have relatively large error 
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Fig. 2. Deuteron analyzing powers \n d -\- d elastic scattering at 
Efi = 10 MeV. Results obtained with potentials INOY04 (solid 
curves), CD Bonn + A (dashed-dotted curves), and CD Bonn (dot¬ 
ted curves) are compared with data from Refs. [38] (•) and [39] (A). 
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Fig. 3. Differential cross section of ^H(d, p)®H (left) and n)®He 

(right) transfer reactions at 10.0 and 12.3 MeV deuteron energy. 
Curves as in Fig. 2. The data are from Refs. [40] (•), [37] (x), [41] 
(□), and [42] (A). 


bars, especially in the iTn case where two data sets [38,39] 
are available. The symmetric analyzing powers T 20 and T 22 
are most sensitive to the employed force model. However, 
the results obtained with different potentials do not corre¬ 
late with the properties of 2N and 3N bound states such as 
binding energies or deuteron U-state probability. The data 
are reproduced best by the predictions using the CD Bonn 
-h A potential. The A-isobar effect is especially significant 
and beneficial for T 20 . The antisymmetric analyzing pow¬ 
ers iTii and T 21 show less sensitivity to the employed po¬ 
tential. At least to some extent this is due to kinematic 
reasons, since iTn are r 2 i vanish exactly at 0c.m. = 90° 
where da/dfl has its minimum and the model dependence 
of symmetric analyzing powers reaches its maximum. 



(Tp 



CD Bonn 

83.4 

87.1 

110 

CD Bonn + A 

79.8 

83.9 

113 

INOY04 

77.6 

81.3 

112 


Table 1 

Predicted total cross sections for ^H(d,p)^H and ^H(d, n)^He trans¬ 
fer reactions and for breakup, labeled Up, (Jn, and cr^, respectively, 
all in millibarns, at 10 MeV deuteron energy. 


4. Transfer reactions 

We calculate the transfer reactions ^H((i,p)^H and 
^H(d, n)^He using the same nuclear interaction models, 
i.e., INOY04, CD Bonn -|- A, and CD Bonn. In Fig. 3 we 
present results at E 4 = 10.0 and 12.3 MeV for the differ¬ 
ential cross section da/dfl as function of the nucleon scat¬ 
tering angle 0c.m. in the c.m. frame. For both p)^H 
and ^H(d, n)^He reactions, da/dVi has very similar shape 
but is slightly higher for the latter. The differential cross 
section is symmetric with respect to 0c.m. = 90° but has 
a more complicated angular and energy dependence than 
in the case of elastic scattering. At both considered ener¬ 
gies there are forward and backward peaks as well as local 
minima around 0c.m. = 45° and 135°. At Ed = 10.0 MeV 
there is just a local maximum located at 0cnL = 90° which 
evolves into a shallow local minimum as the energy in¬ 
creases to 12.3 MeV; meanwhile two local maxima appear 
around 0c.m. = 70° and 110°. Such a behavior is seen also 
in the experimental data [37,40-42]. The overall agreement 
between theoretical results and the data is good, except at 
intermediate angles where da/dU, is small and the data are 
slightly underpredicted by the INOY04 results, roughly 
by 6% at 0c.m. = 90°. The sensitivity to the employed 
potential models is studied at Ed = 10.0 MeV and is vis¬ 
ible around the extrema of da/dVl. While at forward and 
backward peaks and central maximum the predictions 
roughly scale with the trinucleon BE as already observed 
in previous calculations below breakup threshold [24,43], 
at the minima CD Bonn and CD Bonn -|- A results are 
indistinguishable. This may indicate a possible sensitivity 
to two-nucleon isospin singlet partial waves since there CD 
Bonn and CD Bonn -F A potentials are identical. 

The resulting total cross sections ap and Un sX Ed = IQ 
MeV are collected in Table 1 for all three potentials. The 
total transfer cross sections scale quite well with the trin¬ 
ucleon BE. In addition we present also the total breakup 
cross section Uf,, including both three- and four-cluster 
channels. It is calculated using the optical theorem with 
finite screening radius R before subtraction and renormal¬ 
ization of the elastic scattering amplitude because transfer 
and breakup operators are short-ranged and the respec¬ 
tive total cross sections are unchanged by renormalization 
phases [44]. Already at Ed = 10 MeV ab exceeds ap and 
an, indicating the importance of breakup in d + d colli¬ 
sions. Note that in n -F ^He scattering breakup becomes 
the dominant inelastic channel only above the neutron 
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Fig. 4. Deuteron analyzing powers and outgoing nucleon polarization 
of ^H{d,p)^H (left) and ^H(d, n)^He (right) transfer reactions at 
10 MeV deuteron energy. Curves as in Fig. 2. The data are from 
Refs. [40] (+), [45] (□), [46] (x), [47] {■), [48] (o), [49] (.), and [23] 

(A). 

laboratory energy of 23 MeV which roughly corresponds 
to Ed = 28 MeV. 

Next we consider single-polarization spin observables in 
and n)^He reactions at Ed = 10 MeV. In 
Fig. 4 we show vector analyzing power iTn, tensor analyz¬ 
ing powers T 20 , T 21 , and T 22 , and outgoing nucleon polar¬ 
ization Py calculated using INOY04, CD Bonn + A, and 
CD Bonn potentials. All deuteron analyzing powers ex¬ 
hibit a complex angular dependence with several local min¬ 
ima and maxima. They show no symmetry with respect to 
©cm. = 90°, in contrast to Py which is antisymmetric. The 
differences between and ^H(d, n)^He observ¬ 

ables are quite small but visible, e.g., around second max¬ 
imum of T 2 Q. For most observables several measurements 
exist [40,45-49,23] that are at variance in particular cases, 
especially for iTn. Nevertheless, the overall description of 
the experimental data by our calculations is successful with 
only few small or at most moderate disagreements, mostly 
in the vector observables: the minima of iTn are underpre- 
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Fig. 5. Deuteron-to-neutron polarization transfer coefficients of 
^H(d, n)®He reaction at 10 MeV deuteron energy. Curves as in Fig. 2. 
The data are from Ref. [23]. 

dieted while the positive peak of Py is shifted to larger an¬ 
gles. The sensitivity to the used interaction model is quite 
small; it is most visible for iTn. It is not a simple scaling 
with trinucleon BE since often CD Bonn and CD Bonn -|- 
A predictions stay quite close together with INOY04 be¬ 
ing further away. This may again indicate the dominance 
of NN isospin singlet partial waves. It may even indicate 
partial correlation of the observables with the deuteron In¬ 
state probability Pd, but more detailed studies are needed 
to confirm or reject this speculation. 

Finally we show the results for double-polarization ob¬ 
servables. Deuteron-to-neutron polarization transfer coeffi¬ 
cients K^', K^', Kf ,K^', K^', and K^'y in the n)^He 

reaction at Ed = 10 MeV have been measured in Ref. [23]. 
This data and our predictions based on three nuclear in¬ 
teraction models are compared in Fig. 5. The polarization 
transfer coefficients exhibit a very complex angular depen¬ 
dence having up to six local extrema. Given such a compli¬ 
cated behavior of the observables the found agreement be¬ 
tween theory and experiment is impressive. There are only 
small discrepancies such as a slight underestimation of 
at intermediate angles. The model dependence of the polar¬ 
ization transfer coefficients is quite weak and qualitatively 
the same as already seen for single-polarization observables. 

5. Conclusions 

We perform calculations for elastic and transfer reac¬ 
tions initiated by deuteron-deuteron collisions above four- 































nucleon breakup threshold. This process mimics the scat¬ 
tering of two halo nuclei. As dynamic input we use sev¬ 
eral realistic two-nucleon potentials and include the proton- 
proton Coulomb force via the screening and renormaliza¬ 
tion method. Exact four-particle scattering equations in 
the integral form for symmetrized transition operators are 
solved in the momentum-space framework where the pres¬ 
ence of open breakup channels leads to a kernel with a 
highly nontrivial singularity structure. The complex en¬ 
ergy method with special integration weights is successfully 
applied to deal with this difficulty. Compared to previous 
calculations of nucleon-trinucleon scattering, the relatively 
weak binding of deuteron and its large spatial size lead 
to additional complications such as a slower partial-wave 
and complex-energy convergence. Nevertheless, we obtain 
fully converged results for d -|- d elastic scattering as well 
as for ^H(d,p)^H and ^H(d, n)^He transfer reactions. For 
these reactions at 10 MeV deuteron energy we calculate the 
differential cross section and all deuteron analyzing pow¬ 
ers; the former observable is predicted also at Ed = 12.3 
MeV. Furthermore, for transfer reactions we calculate also 
the outgoing nucleon polarization and, in the ^H(d, n)^He 
case, deuteron-to-neutron polarization transfer coefficients. 
The overall description of the experimental data is good, 
even for the most complicated double-polarization observ¬ 
ables. The comparison of predictions based 1NOY04, CD 
Bonn -|- A, and CD Bonn potential models may indicate 
the dominance of NN isospin singlet partial waves for most 
spin observables in transfer reactions, but not in the case 
of elastic scattering. We also predicted the total breakup 
cross section and demonstrated the increased importance 
of breakup channels in d + d reactions. 

Together with the previous achievements in the nucleon- 
trinucleon scattering, the present work demonstrates that 
numerically exact calculations of all two-cluster reactions 
in the four-nucleon system are now possible in a fully con¬ 
verged way using realistic nuclear interactions and includ¬ 
ing the pp Coulomb repulsion. 
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